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Entanglement of two parts of a quantum system is a nonlo- 
cal property unaffected by local manipulations of these parts. 
It is described by quantities invariant under local unitary 
transformations. Here we present, for a system of two qubits, 
a set of invariants which provides a complete description of 
nonlocal properties. The set contains 18 real polynomials of 
the entries of the density matrix. We prove that one of two 
mixed states can be transformed into the other by single- 
bit operations if and only if these states have equal values 
of all 18 invariants. Corresponding local operations can be 
found efficiently. Without any of these 18 invariants the set 
is incomplete. 

Similarly, nonlocal, entangling properties of two-qubit uni- 
tary gates are invariant under single-bit operations. We 
present a complete set of 3 real polynomial invariants of uni- 
tary gates. Our results are useful for optimization of quantum 
computations since they provide an effective tool to verify if 
and how a given two-qubit operation can be performed us- 
ing exactly one elementary two-qubit gate, implemented by a 
basic physical manipulation (and arbitrarily many single-bit 
gates) . 

Introduction. Nonlocality is an important ingredient 
in quantum information processing, e.g. in quantum 
computation and quantum communication. Nonlocal 
correlations in quantum systems reflect entanglement be- 
tween its parts. Genuine nonlocal properties should be 
described in a form invariant under local unitary oper- 
ations. In this paper we discuss such locally invariant 
properties of (i) unitary transformations and (ii) mixed 
states of a two-qubit system. 

Two unitary transformations (logic gates), M and L, 
are called locally equivalent if they differ only by local 
operations: L = UiMU 2 , where Ui,U 2 £ SU{2) m are 
combinations of single-bit gates on two qubits Q. A 
property of a two-qubit operation can be considered non- 
local only if it has the same value for locally equivalent 
gates. We present a complete set of local invariants of 
a two-qubit gate: two gates are equivalent if and only 
if they have equal values of all these invariants. The 
set contains three real polynomials of the entries of the 
gate's matrix. This set is minimal: the group SU (4) 
of two-qubit gates is 15-dimcnsional and local operations 
eliminate 2 dim[Sf7 (2)® 2 ] = 12 degrees of freedom; hence 
any set should contain at least 15—12 = 3 invariants. 

This result can be used to optimize quantum compu- 
tations. Quantum algorithms are built out of elementary 
quantum logic gates. Any many-qubit quantum logic cir- 
cuit can be constructed out of single-bit and two-bit op- 
erations |2|-||. The ability to perform 1-bit and 2-bit 



operations is a requirement to any physical realization 
of quantum computers. Barenco et al. [0 showed that 
the controlled-not (CNOT) gate together with single-bit 
gates is sufficient for quantum computations. Further- 
more, it is easy to prove that any two-qubit gate M forms 
a universal set with single-bit gates, if M itself is not 
a combination of single-bit operations and the SWAP- 
gate, which interchanges the states of two qubits. The 
efficiency of such a universal set, that is the number of 
operations from the set needed to build a certain quan- 
tum algorithm, depends on M. 

For a particular realization of quantum computers, a 
certain two-qubit operation M (or a set exp[ffl.t] of op- 
erations generated by a certain hamiltonian) is usually 
considered elementary. It can be performed by a basic 
physical manipulation (switching of one parameter or ap- 
plication of a pulse). Then the question of optimization 
arises: what is the most economical way to perform a 
particular computation, i.e. what is the minimal num- 
ber of elementary steps? In many situations two-qubit 
gates are more costly than single-bit gates (e.g., they can 
take a longer time, involve complicated manipulations 
or stronger additional decoherence), and then only the 
number of two-bit gates counts. 

The simplest and important version of this question is 
how a given two-bit gate L can be performed using the 
minimal number of elementary two-bit gates M. In par- 
ticular, when is it sufficient to employ M only once 
This is the case when the two gates are locally equivalent, 
and computation of invariants gives an effective tool to 
verify this. Moreover, if M and L are equivalent, a proce- 
dure presented below allows to find efficiently single-bit 
gates Ui, U2, which in combination with M produce L. 
If one elementary two-bit step is not sufficient, one can 
ask how many are needed. Counting of dimensions sug- 
gests that two steps always suffice. However, this is not 
true for some gates M (bad entanglers). 

A related problem is that of local invariants of quan- 
tum states. A mixed state is described by a density ma- 
trix p. Two states are called locally equivalent if one 
can be transformed into the other by local operations: 
p — > W pU, where U is a local gate. Apparently, the co- 
efficients of the characteristic polynomials of p and of the 
reduced density matrices of two qubits are locally invari- 
ant. The method developed in Refs. in principle, 
allows to compute all invariants. For a two-qubit sys- 
tem counting of dimensions shows that there should be 
9 functionally independent invariants of density matrices 
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with unit trace, but additional invariants may be needed 
to resolve a remaining finite number of states. A set of 
20 invariants was presented in Ref. B). However it was 
not clear if this set was complete, i.e. two states with the 
same invariants are always locally equivalent. Here we 
present a complete set of 18 polynomial invariants. We 
prove that two states are locally equivalent if and only if 
all 18 invariants have equal values in these states. Hence, 
any nonlocal characteristic of entanglement is a function 
of these invariants ||. We also show that no subset of 
this set is complete. 

To demonstrate applications of our results, we discuss 
in the last section which 2-bit operations can be con- 
structed out of only one elementary 2-bit gate for Joseph- 
son charge qubits |ic],[llj or for qubits based on spin de- 
grees of freedom in quantum dots 

Single-bit gates as orthogonal matrices. The following 
result is used below to classify two-qubit gates. 

Theorem 1. Single-qubit gates with unit determinant 
are represented by real orthogonal matrices in the Bell 
basis ^(|00) + |11)), ^(|01) + |10>), ^(|01) - |10», 

^(|00) -|11)). 

The transformation of a matrix M from the standard 
basis of states 1 00) , 1 01) , 1 10) , |11) into the Bell basis is 
described as M —> Mb = Q^MQ, where 
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Proof. We introduce a measure of entanglement in a 
pure 2-qubit state ip a /3, & quadratic form EntV'j which 
in the standard basis is defined as Ent ip = det xp — 
V'ooV'ii — V'oiV'io- This quantity is locally invariant. 
Indeed, a single-bit gate Wi ® Wi transforms ip into 
\V1ipW2 , preserving the determinant. 

Under a unitary operation V the matrix of this form 
transforms as Ent — > U T Ent V. In the Bell basis it 
is proportional to the identity matrix. Since local gates 
preserve this form, they are given by orthogonal matrices 
in this basis. As unitary and orthogonal they are also 
real. Thus local operations form a subgroup of the group 
50(4, R) of real orthogonal matrices. This subgroup is 
6-dimensional, and hence coincides with the group. □ 

Classification of two-qubit gates. Our result for uni- 
tary gates is expressed in terms of Mb ■ 

Theorem 2. The complete set of local invariants of a 
two-qubit gate Af , with det M — 1, is given by the set of 
eigenvalues of the matrix M b Mb- 

In other words, two 2-bit gates with unit determinants, 
M and L, arc equivalent up to single-bit operations iff 
the spectra of M b Mb and L^Lb coincide. Since m = 
M^Mb is unitary, its eigenvalues have absolute value 
1 and are bound by detm = 1. For such matrices the 



spectrum is completely described by a complex number 
tr to and a real number tr 2 777 — tr m 2 . 

The matrix m is unitary and symmetric. The following 
statement is used in the proof of Theorem 2: 

Lemma. Any unitary symmetric matrix m has a real 
orthogonal eigenbasis. 

Proof. Any eigenbasis of 777 can be converted into a real 
orthogonal one, as seen from the following observation: If 
v is an eigenvector of m with eigenvalue A then v* is also 
an eigenvector with the same eigenvalue (conjugation of 
777*v = 777 _1 v = A _1 v = A*v gives this result). Hence, 
Rev and Imv are also eigenvectors. □ 

Proof of Theorem 2. In the Bell basis local opera- 
tions transform a unitary gate Mb into O1M5O2, where 
Oi,0 2 € 50(4, R) are orthogonal matrices. Therefore 
777 = M^Mb is transformed to OJ777O2. Obviously, the 
spectrum of 777 is invariant under this transformation. 

To prove completeness of the set of invariants, we no- 
tice that the lemma above implies that m can be di- 
agonalized by an orthogonal rotation Om, i.e. m = 
O^cImOm where oIm is a diagonal matrix. Suppose that 
another gate L is given, and I = L^Lb has the same 
spectrum as m. Then the entries of oIm and d^ are re- 
lated by a permutation. Hence cLm = P T dr 1 P, where P is 
an orthogonal matrix, which permutes the basis vectors. 
Using the relation of m to c?m and of I to dr,, we conclude 
that I = Q T mO where O € 50(4, R). 



Single-bit operations O and O' = LbO M b 
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form one gate into the other: Lb = O'MbO. The 
gate O' is a single-bit operation since it is real and or- 
thogonal. Indeed, 0' T 0' = {M B X ) T OL T ? L B T M B X = 
(M B 1 ) T OlO T M B 1 = (Mg^mMg 1 = 1. On the other 
hand, O' is unitary as a product of unitary matrices. This 
implicates its reality. □ 

So far we have discussed equivalence up to single-bit 
transformations with unit determinant. However, physi- 
cally unitary gates are defined only up to an overall phase 
factor. The condition det M — 1 fixes this phase factor 
but not completely: multiplication by ±7 preserves the 
determinant. With this in mind, we describe a procedure 
to verify if two two-qubit unitary gates with arbitrary de- 
terminants are equivalent up to local transformations and 
an overall phase factor: 

We calculate m = M b Mb for each of them and com- 
pare the pairs [tr 2 m det Aft ; trm 2 det Aft]. If they co- 
incide, the gates are equivalent and the proof of Theorem 
2 allows to express explicitly one gate via the other and 
single-bit gates, O and O 1 . 

Classification of two-qubit states. In this section we 
discuss equivalence and invariants of two-qubit states up 
to local operations. Let us express a 2-bit density matrix 
in terms of Pauli matrices acting on the first and the 



second qubit: p = -A. 



tpa 2 + fajaja? If the 



qubits are considered as spin-1/2 particles, then s and p 
are their average spins in the state p, while (3 is the spin- 
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spin correlator: /3jj = (S^iS 2 ). Any single-bit operation 
is represented by two corresponding 3x3 orthogonal real 
matrices of 'spin rotations', 0,P G 50(3, R). Such an 
operation, O ® P, transforms p according to the rules: 
s -► Os, p -> Pp, and /3 -> 0/3P T . We find a set of 
invariants which completely characterize p up to local 
gates. 

The density matrix is specified by 15 real parameters, 
while local gates form a 6-dimensional group. Thus we 
expect 15 — 6 = 9 functionally independent invariants 
(h-Ig in Table |). However, these invariants fix a state 
only up to a finite symmetry group, and additional in- 
variants are needed ||. In Table I we present a set of 18 
polynomial invariants and prove that the set is complete. 
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TABLE I. The complete set of invariants of a two-qubit 
state. Here (a, b, c) stands for the triple scalar product 
a - (b x c), and eijk is the Levi-Cevita symbol; s/3 is a 3- vector 
with components Sj/3jj etc. The last two columns present, 
for each invariant, s, p in a situation when only this invariant 
distinguishes between two nonequivalent states. In these ex- 
amples we assume a nondegenerate /3, with &3 = for Ji2-i8- 

Theorem 3. Two states are locally equivalent exactly 
when the invariants h~hs have equal values for these 
states. (Only signs of iio, In, lis-ig are needed.) 

None of the invariants can be removed from the set 
without affecting completeness, as demonstrated by ex- 
amples in the table. 

The proof below gives an explicit procedure to find 
single-bit gates which transform one of two equivalent 
states into the other. 

Proof. It is clear that all I in the table are invariant 
under independent orthogonal rotations O, P, i.e. under 
single-bit gates. To prove that they form a complete set, 
we show that for given values of the invariants one can 
by local operations transform any density matrix with 
these invariants to a specific form. In the course of the 
proof we fix more and more details of the density matrix 
by applying local gates (this preserves the invariants). 

The first step is to diagonalize the matrix j3, which can 
be achieved by proper rotations O, P (singular value de- 



composition). The invariants I1—I3 determine the diago- 
nal entries of (3 up to a simultaneous sign change for any 
two of them. Using single-bit operations K 1 ® 1 (where 
R l is the 7r-rotation about the axis i = 1,2, 3) we can fix 
these signs: all three eigenvalues, 61,62, 63, can be made 
nonnegative, if h = det $ > 0, or negative, if det (3 < 0. 
Further transformations, with O = P representing per- 
mutations of basis vectors, place them in any needed or- 
der. From now on we consider only states with a fixed 
diagonal (3. Hence, only such single-bit gates, 0®P, are 
allowed which preserve j3. The group of such operations 
depends on b\,b 2 ,b 3 . We examine all possibilities below, 
showing that the invariants h—xs fix the state completely. 

(A) (3 is nondegenerate: all 6, are different. Then the 
only /3-preserving local gates are r % = R l £g> i?\ 

The invariants hs and I7-9 set absolute values of six 
components Sj, pi, but not their signs. The signs are 
bound by other invariants. In particular, ho and In fix 
the values of S1S2S3 andpip 2 p 3 . Furthermore, /12-14 give 
three linear constraints on three quantities SiPi. When 
(3 is not degenerate, one can solve them for SiPi. Let us 
consider several cases: 

i) s has at least two nonzero components, say s\,S2- 
These two can be made positive by single-bit gates r 1 , 
r 2 . After that, the signs of Pi=i,2 are fixed by values of 
SiPi, while iio fixes the sign of S3. The sign of p 3 can be 
determined from s 3 p 3 , if S3 7^ 0; if S3 = then p$ is fixed 
by hs = P3S 1 s 2 b 3 (bl - b\) or I 17 = p 3 s 1 s 2 b 1 b 2 {bl - b\). 

If p has at least two nonzero components, a similar 
argument applies, with /11,16,1s instead of Iio, 15, 17- If 
both s and p have at most one nonzero component, s, 
and pj, then either ii) i = j, and the signs are specified 
by SiPi up to /3-preserving gates r k ; or iii) i ^ j, and one 
can use r l ,ri to make both components nonnegative. 

(B) (3 has two equal nonzero eigenvalues: 63 7^ b\ = 
62 7^ 0. We define horizontal components = (s±, s 2 , 0), 
p^ = (pi,p 2 ,0). Then /3-preserving operations are gen- 
erated by simultaneous, coinciding rotations of Sj^ and 
Px [i.e. Sx — > Osj_, pj^ — > Op_L, where O S 501^(2) is 
a 2 D- rot at ion] , as well as r % . 

The invariants 74_g fix Sj_, s|, p^ and p\. To specify s 
and p completely the angle between sj^ and p^, as well 
as the signs of S3 and p 3 should be determined. These are 
bound by the remaining invariants. In particular, /12-14 
fix s 3 p 3 and s_i_pj_. The latter sets the angle between sj_ 
and p^ up to a sign. 

There are two possibilities: i) S3 = p 3 = 0. The states 
with opposite angles are related by r 1 and hence equiva- 
lent, ii) s 3 7^ (the case p 3 7^ is analogous). Applying 
r , if needed, we can assume that S3 is positive. Then, 
p 3 is specified by the value of s 3 p 3 . Apart from that, hb 
sets (sj_ x PJJ3S3, and hence the sign of the cross product 
s±_ x p^. This fixes the density matrix completely. 

(C) b\ = b 2 = 0, 63 7^ 0. In this case /3-preserving 
operations are independent rotations of sj_ and p^ [which 
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form SOi^i^)® 2 ] and r\ The invariants I±-g fix Sj_, s\, 
p^, and p|, while /12 (or I 13 ) sets S3P3. It is easy to see 
that they specify the state completely. 

(D) 61 = 62 = ^3 7^ 0. All transformations O <8> O, 
where O € SO(3), preserve j3. The invariants J4, J7 and 
J12 fix s 2 , p 2 and sp, and this information is sufficient to 
determine s and p up to a rotation. 

(E) /3 = 0. In this case all local transformations pre- 
serve j3\ hence, s and p can rotate independently. The 
only nonzero invariants, I4 and J7, fix s 2 and p 2 . □ 

Discussion. In this section we demonstrate applica- 
tions of our results. We calculate the invariants for sev- 
eral two-qubit gates to find out which of them are locally 
equivalent. These gates include CNOT, SWAP and its 
square root, as well as several gates exp(iH.t) generated 
by hamiltonians ja 1 ^ 2 , j<j]_a^_ [here a± = (a x ;a y )] and 
ifiyCr 2 (cf. Refs. [ pr| , ^2| J^1) after evolution during time t. 



%/SWAP are perfect entanglers. The Heisenberg hamilto- 
nian can produce only two perfect entanglers (VSWAP or 
its inverse), while (J v <J y - only CNOT. At the same time, 
the interaction a±a± produces a set of perfect entanglers 
if the system evolves for time t with cost < 0. 

To conclude, we have presented complete sets of lo- 
cal polynomial invariants of two-qubit gates (3 real in- 
variants) and two-qubit mixed states (18 invariants) and 
demonstrated how these results can be used to optimize 
quantum logic circuits and to study entangling properties 
of unitary operations. 

I am grateful to G. Burkard, D.P. DiVincenzo, 
M. Grassl, G. Schdn, and A. Shnirman for fruitful dis- 
cussions. 
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TABLE II. Invariants of two-qubit gates, 
Gi = tr 2 m/16 det M and G 2 = (tr 2 m - tr m 2 )/4 det M. The 
latter is always a real number. Numerical prefactors are cho- 
sen to simplify expressions. 



Analyzing the invariants we see that to achieve CNOT 
one needs to perform a two-qubit gate at least twice if 
the latter is triggered by the Heisenberg hamiltonian aa. 
At the same time SWAP and -\/SWAP can be performed 
with one elementary two-bit gate ||. The interaction 
<J y <Jy allows to perform CNOT in one step, while a±a± 
requires at least two steps for all three gates. 

In Josephson qubits jlf| elementary two-bit gates are 
generated by the interaction H = — ^Ej(a^ + a 2 ) + 
(Ej/E^fry&y. Investigation of the invariants shows that 
CNOT can be performed if Ej is tuned to aE^ for a fi- 
nite time t — aTr(2n + 1)/4El, where n is an integer and 
a satisfies a 2 cos [71- (n + + a~ 2 ] = — I . 

For creation of entanglement between qubits a useful 
property of a two-qubit gate is its ability to produce a 
maximally entangled state (with | Ent?/>| = 1/2) from an 
unentangled one |J. This property is locally invariant 
and one can show that a gate M is a perfect entangler 
exactly when the convex hull of the eigenvalues of the 
corresponding matrix m contains zero. In terms of the 
invariants, introduced in Table 0, this condition reads: 
sin 2 7 < 4|Gi| < 1 and cos7(cos7 — G2) > 0, where 
d = |Gi|e 17 . Among the gates in the table CNOT and 
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